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Abstract 



We present a fully dynamic algorithm for maintaining approximate maximum weight 
matching in general weighted graphs. The algorithm maintains a matching M whose size 
is at least \M* where M* is the weight of the maximum weight matching. The algorithm 
achieves an expected amortized O(lognlogC) update time, where C is the ratio between 
the maximum and the minimum weights of edges in the given graph. Using a simple 
randomized scaling technique, we are able to obtain a matching whose expected size is at 
least \M*. 

1 Introduction 

Given an undirected graph G = (V, E), where V represents the set of vertices and E represents 
the set of edges in the graph, let n = \V\ and m = \E\. In a weighted graph, we associate a 
weight function w : E — > M where Ve G E, w(e) represents the weight of the edge. The weight 
function for a set of edges M C E is represented by w(M) and defined as Yle&M w ( e )- 

A matching in a graph is a set of independent edges in the graph. Thus, a subset M of 
E is called a matching if no vertex of the graph is incident on more than one edge in M. In 
an unweighted graph, we refer maximum matching as maximum cardinality matching(MCM) 
while in a weighted graph we refer to it as a maximum weight matching(MWM). A matching 
is called o-MWM(o-MCM) if it is at least ± factor of MWM(MCM), and is represented by 
a-MWM (a-MCM) . 

A dynamic graph algorithm aims to answer queries regarding some property of the graph 
which is undergoing a series of changes due to insertion or deletion of edges (V is assumed to 
be fixed). The updates can be classified naturally into two types - insertions and deletions. 
Algorithms which handle only insertions are called incremental and the ones that can deal 
only with deletions are called decremental. An algorithm that can handle both insertions and 
deletions of edges is called a fully dynamic algorithm. In this paper, we address the problem 
of maintaining approximate weighted matching in a fully dynamic graph. 



Previous Results 

The fastest known algorithm for finding MCM in general graphs is by Micali and Vazirani[7J 
that runs in Oiyn^fn) time. Their algorithm can be used to compute a matching having size 
(1 — e) times the size of maximum matching in 0{m/e) time. Mucha and Sankowski[8j gave 
an algorithm that computes MCM in 0(n?), where u < 2.376 is the exponent of n in the 
fastest known matrix multiplication algorithm. Relatively few algorithms are known that 
maintain matching in dynamic graphs. The first result was presented by Ivkovic and Lloyd[5j 
with amortized update time 0(n + m) ' 7072 . Onak and Rubinfeld|9] achieved polylogarith- 
mic amortized expected update time and maintained an a-approximate MCM where a was 
shown to be some constant but not explicitly calculated. Baswana, Gupta and Sen[T] pre- 
sented a fully dynamic randomized algorithm for maintaining maxima^] matching in expected 
amortized O(logn) update time. 

For MWM, Gabow[4j gave an 0(mn + ra 2 logra) time algorithm. Preis[10] gave a 0(m) 
time algorithm for finding 2- MWM. Drake and Hougardy[2] gave a simpler algorithm for the 
same problem. Vinkemeier and Hougardy[TT] presented an algorithm to get a matching which 
is (2/3 — e) times the size of MWM in 0(m/e) time. Duan, Pettie and Su[3] presented an 
algorithm to find a matching which is (1 — e) times the size of MWM in 0(me _1 log e _1 ) time. 
To the best of our knowledge, there have been no significant results for maintaining MWM 
or approximate MWM in dynamic graphs. 

Preliminaries 

Let M be the set of matched edges in the graph. A vertex is called free if it is not incident 
on any edge in M. A vertex which is not free is called matched. Similarly, an edge is called 
matched if it is in M and is called free otherwise. Suppose (u, v) is a matched edge, then u is 
said to be the mate of v and vice versa. A matching M is called maximal if no edge can be 
added to the matching without violating the degree bound of one for a matched vertex. An 
alternating path is defined as a path in which edges are alternately matched and free, while 
an augmenting path is an alternating path which begins and ends with free vertices. 

Our Results 

We present a fully dynamic algorithm that achieves expected amortized O(lognlogC) update 
time for maintaining 8-MWM. Subsequently, we improve this approximation ratio to 6 using 
a simple randomized scaling technique. Here C is the ratio of the highest weight edge to the 
smallest weight edge in the given graph. It uses the result of Baswana, Gupta and Sen [I] as 
a subroutine. The algorithm uses some ideas from an approximate (distributed) algorithm 
given by Lotker, Patt-Shamir, and Rosen[6j. They give a distributed algorithm having running 
time 0(1) per update and which mainains a 27-MWM in a dynamic weighted graph. The 
first subsection describes the basic ideas and the intuition behind the algorithm. Section 
describes our algorithm. 

Overview of our approach 

The basic idea of the algorithm is to try and reduce the problem of finding an approximate 
maximum weight matching to the problem of finding approximate matching in unweighted 

1 & 2-MCM 



2 



graphs. Observe that finding MCM is a special case of finding MWM with all edges having 
the same weight. So, a maximal matching is 2-MWM when all edges have the same weight. 
But this observation does not immediately extend to general the weights. First, let us con- 
sider the case when the weights are within a range, say, [a 1 , a t+1 ), where a > 1 is a constant. 
In such a graph the maximal matching gives a 2a approximation of the maximum weight 
matching. This follows directly from the fact that maximal matching in an undirected graph 
gives 2 approximation of maximum matching. So, maximal matching can be used as an ap- 
proximation for MWM in a graphs where the weights of the edges are bounded. To generalize 
this, we partition the edges of the graph into levels according to their weight. Edges at level 

1 have weights in the range [a 1 , a l+1 ) and the set of edges at level i is represented by Ei, viz., 
Ve G Ei,w{e) G [a\a i+l ). 

Observe that in this scheme of partitioning, any edge is present only at one level while a 
vertex may have edges incident on it from various levels. The subgraph at level i is represented 
by Gi = (V,Ei) and a maximal matching Mi is maintained for G{. The maximal matching 
at each level provides an approximation for the maximum matching at that level. To obtain 
an overall bound on MWM, the approximations obtained by the maximal matchings from 
different levels have to be combined in some way. The need for this combining arises from 
our earlier observation that even if only the edges in UjMj are considered, a vertex may have 
multiple edges incident on it. So, some of the edges have to be removed from UjMj to find a 
valid matching of G. Let H = (V, (J Mj) be the subgraph of G having only those edges which 
are part of the maximal matching at some level. In this subgraph, a matching is maintained 
which is shown to be 8-MWM. 

2 Fully Dynamic 8-MWM 

We maintain a partition of edges according to their levels. A maximal matching Mj is main- 
tained at each level. We use the fully dynamic algorithm in [1] to maintain Mj's. Their 
algorithm will cause edges to become matched from free and vice-versa leading to addition 
or deletion of edges to the subgraph %. They show that there is an algorithm which can 
maintain maximal matching in a dynamic graph in O(logn) amortized time with high prob- 
ability. This means that the amortized number of matched edges added or removed from the 
subgraph % is O(logn) at each update step. We maintain the output matching A4 in this 
subgraph % taking advantage of the hierarchical structure of T~L. Since T~L is formed by the 
union of matchings at various levels, a vertex can have at most one neighbor at each level. 
The matching is maintained such that for every edge of T~L which is not in A4 there is an edge 
adjacent to it at a higher level which is in Ai. For an edge e, let Level(e) denote its level. 
Formally the invariant is, 

Ve G E(7i), either e G M. or e is adjacent to an edge e' G M. such that Level (e') > 
Level{e). 

Data Structures Used 

The algorithm requires information to be maintained about the graph and its vertices and 
edges. We explicitly show how the following are maintained - 

• Mi - A maximal matching at the level I. 
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• Free(v) - Whether the vertex v is free in the matching M.. 

• Mate(v) - The mate of v, if it is not free. 

• Level((u,v)) or Level{e) - The level at which the edge e or the edge (u,v) is according 
to the condition that Ve G Gi,w(e) G [a l ,a* +1 ). 

• OccupiedLevels - The set of levels where there is at least one edge from %. 

• L max - The highest occupied level. 

• L mm - The lowest occupied level. 

• N(v, I) - It maintains the neighbor of v at the level I if there is one, null otherwise. 

• A4 - The matching being maintained. 

The next subsection describes a static algorithm for computing M. in the graph %. The 
static algorithm gives way to a dynamic algorithm described in the subsection following it. 

Static Algorithm to obtain Ai from % 
Procedure 2.1: StaticCombine () 



The static algorithm divides the edges from the graph G into levels and a maximal match- 
ing Mi is obtained for each of the levels. Using these maximal matchings we get the graph %. 
Thereafter the level numbers L max and L mm are computed and the procedure StaticCombine 
is used. 

The procedure StaticCombine starts by picking all the edges in H at the highest level and 
adds them to the matching M. For every edge (u, v) added to the matching M, all the edges 
in the graph % incident on u and v have to be removed from the graph. The same process is 
repeated for the next lower level. Note that every edge in % is either in the matching M or 
is deleted due to the addition of a neighboring edge in M. and thus maintains the invariant. 
Also, the matching M. is maximal in H because of the way it is being computed. 



1 M = <p; 

2 for i = L max to L min do 

3 M=MUMi; 

4 for (u, v) G Mi do 

5 for j = i - 1 to L min do 

6 for (x, y) G Mj do 

7 if u = x or u = y or v = x or v = y then 
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Dynamic Algorithm to maintain Ai 

Before any addition or deletion of edge, we assume that there was a matching M which 
satisfied the invariant. We begin the description of the algorithm in terms of addition and 
deletion of edges to the graph ri. The addition and deletion of the edges in ri is caused due 
to addition/deletion of an edge in G. We describe some of the frequently used procedures 
first. Then the procedures for addition and deletion of edges in % and finally the procedures 
for addition and deletion of edges in G are described. 



Procedure 2.2: AddToMatching(w, v) 

1 Free(u) = False; Free(v) = False; 

2 Mate(u) = v; Mate(v) = u; 

3 M = M\J{(u,v)}; 



Procedure 2.3: DelFromMatching(u, v ) 



1 Free(u) = True; Free(v) 

2 M = M\{(u,v)}; 



True; 



The procedure AddToMatching adds an edge to the matching updating the free and mate 
fields accordingly. The procedure DelFromMatching deletes an edge from the matching up- 
dating the mate and the free fields correctly. Both of them execute in O(l) time. 



Procedure 2.4: HandleFree(«, lev ) 



l for I from lev to L mm do 



2 
3 
4 
5 
6 

7 
8 
9 
10 
11 
12 



v = N(u,l); 
if v is not null then 
if v is free then 

AddToMatching (u, v); 
return; 

else if Level ((v, Mate(v))) < I then 
v > = Mate{v); 
DelFromMatching (v,v'); 
AddToMatching (u,v); 
HandleFree (v' ', Level((v,v'))); 
return; 



The procedure HandleFree takes as an input a vertex u which has become free and a level 
number lev from where it has to start looking for a mate. The procedure assumes that u 
cannot find a mate at any level above lev. The procedure checks the neighbors of the vertex 
u at all levels below lev falling step by step. If it finds a free neighbor the corresponding edge 
is added to the matching and thereafter nothing is done. Otherwise if a vertex is found which 
has a mate at a lower level then we are violating the invariant as (u, v) does not belong to 
M and is neighboring to an edge in M at a lower level. So, the edge (v, Mate(v)) is removed 



5 



from the matching and the edge (u, v) is added to the matching. This change results in a 
free vertex which is at a lower level which has to be handled. Note that the recursive calls to 
HandleFree are all with lower level numbers. So, the procedure takes 0{L max — L mm ) time. 



Procedure 2.5: AddEdge(u, v) 



1 I = Level((u, v)); 

2 N(u,l) = v; 

3 N(v, I) = u; 

4 if u is free and v is free then 

5 |_ AddToMatching (u,v); 

6 else if u is free and v is not free then 

7 if Level((u, Mate(u))) < I then 

8 v ' = Mate(v); 

9 DelFromMatching (v , v'); 
10 AddToMatching (u, v); 

n HandleFree (v', Level((v, v'))); 

12 else if u is not free and v is free then 

13 if Level((v, Mate(v))) < I then 

14 v! = Mate(u); 

15 DelFromMatching (u, u'); 

16 AddToMatching (u, v); 

17 HandleFree (u' , Level ((u, u'))); 

is else if Level((v, Mate(v))) < I and level((u, Mate(u))) < I then 

19 
20 
21 
22 
23 



u' = Mate(u); v' = Mate(v); 

DelFromMatching (u, u'); DelFromMatching (v,v'); 
AddToMatching (u,v); 
HandleFree {v! , Level((u,u'))); 
HandleFree (v', LevelUv , v'))); 



The procedure AddEdge handles addition of edges to "H. Suppose the edge (u, v) is added 
to %. It is checked whether both u and v are free or not. If so the edge (u,v) is added to 
the matching M. If not then at least one of u and v is matched. It checks for the edge (u, v) 
whether it is adjacent to an edge in M which is at a higher level or not. If (u,v) is adjacent 
to a higher level edge in Ai, then nothing is done else (u, v) is adjacent to lower level edge(s) 
in Ai, thus violating the invariant. The procedure removes the lower level edge(s) from the 
matching and adds the edge (u, v) to the matching. At most 2 vertices can become free due 
to the addition of this edge which are handled using the procedure HandleFree. If v! was the 
previous mate of u, the the edge (u,u') is removed from Ai. Since Ai satisfied the invariant 
before addition of this edge, all the neighboring edges of v! at higher level than Level(u,u') 
are matched to a vertex at higher levels. So v! has to start looking for mates from the level 
of (u,u f ). The procedure makes a constant number of calls to HandleFree and thus runs in 
0(L max - L min ) time. 

The procedure DeleteEdge does nothing if an unmatched edge is deleted. If a matched 
edge (u,v) is deleted at level I, it calls HandleFree for both the end points to restore the 
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Procedure 2.6: DeleteEdge(it, v) 



1 I = Level((u, v)); 

2 N(u, I) = null; 

3 N(v, I) = null; 

4 if (u, v) G M then 

5 DelFromMatching (u, v ); 

6 HandleFree (u, I); 

7 HandleFree (v, I); 



invariant in the matching. HandleFree is called with the level I because our invariant implies 
that all the neighbors of u and v are matched at higher levels. So they cannot find a mate at 
higher levels. This again takes 0(L max — L mm ) time. 



Procedure 2.7: EdgeUpdate(u, v, type) 



1 I = Level((u,v)) = [log a w(u,v)\; 

2 if type is addition and Mi is 4> then 

3 OccupiedLevels = OccupiedLevels\J{l}; 

4 |_ Update L max and L min ; 

5 Update Mi using the algorithm in [1]; 

6 if type is deletion and Mi is 4> then 

7 OccupiedLevels = OccupiedLevels \ {I}; 

8 |_ Update L max and L min ; 

9 Let V be the set of edges deleted from Mi in step 5; 
10 Let A be the set of edges added to M; in step 5; 

n for (x, y) G V do 

12 |^ DeleteEdge (x, y); 

13 for (x, y) G A do 

14 |_ AddEdge (x,y); 



The function EdgeUpdate handles addition and deletion of edge in G. It finds out the level 
of the edge and updates the maximal matching at that level. It updates the OccupiedLevels 
set accordingly. This set is required because the values of L max and L min are to be maintained. 
The algorithm in [1] can be easily augmented to return the set of edges being added or deleted 
from the maximal matching in each update. As discussed before, amortized O(logn) edges 
change their status per update. So, the procedure and the algorithm has an expected amor- 
tized update time of 0(logn • {L max — L mm )). Let e max and e mm represent the edges having 
the maximum and the minimum weight in the graph. Recall that C = w(e max )/w(e mm ). 



L ma X _ L rmn < ^ ^ma^ _ ^ + 1 = Q , bg _±__J_ ) = Q{]ogC) 

So we can claim that 
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Claim 2.1. The expected amortized update time of the algorithm per edge insertion or deletion 
is O(lognlogC). 

In the next section we analyze the algorithm to prove that there is a choice of a for which 
the approximation factor becomes 8. 

2.1 Analysis 

To get a good approximation ratio, we bound the weight of M* with the weight of Ai. We 
now state a few simple observations which help in understanding the analysis. 

Observation 2.1. Since M* is a matching, no two edges of M* can be incident on the same 
vertex. 

Observation 2.2. For any edge e £ M* , there can be at most two edges of M* which are 
adjacent to e, one for each endpoint of e. 

To bound the weight of M* using the weight of Ai, we define a many to one mapping 
(p : E — > E. This mapping maps every edge in M* to an edge in Ai. Using this mapping, we 
find out all the edges which are mapped to an edge e G Ai and bound their weight using the 
weight of e. Let this set be denoted by _1 (e). For an edge e G M*, the mapping is defined 
as: 

1. If e G E{U) and e G Ai then cp(e) = e. 

2. If e G E(T~L) and e ^ Ai then our invariant ensures that e is adjacent to an edge e' such 
that Level(e') > Level(e), then 0(e) = e'. If e is adjacent to two matched edges in A4, 
map e to any one of them. As a rule, if two edges are available for mapping, then we 
will map e to any one of them. 

3. If e ^ E{%), then e is adjacent to an edge e' G E(Ti), if e' G Ai then </>(e) = e'. 

4. If e ^ E(H), then e is adjacent to an edge e' G E(7i), if e' ^ A1 then e' is adjacent to 
an edge e" G .M such that Level(e") > Level(e'), c/>(e) = e". 

Now that we have defined a many to one mapping, we find out the edges of M* which are 
mapped to an edge e G M. An edge which is mapped to e can either be adjacent to e or e 
itself or not adjacent to e. If an edge of M*, which is mapped to e G Ai, is adjacent to e or 
is e itself, then it is called a Directly mapped edge. An edge of M* which is mapped to e G 
M. and is not adjacent to e is called an Indirectly mapped edge. Let 0^ 1 (e) and <i>J l {e) be the 
set of directly mapped and indirectly mapped edges respectively for an edge e G M. Directly 
mapped edges are of type 1, 2 and 3 as described above and indirectly mapped edges are of 
type 4. Since we maintain a maximal matching Ai, this ensures that an edge e G M* is in 
A4 or is adjacent to an edge in A4. This ensures that all the edges in M* are mapped by 4>. 

A directly mapped edge maybe of the first type as described above. If an edge e G M. 
has a directly mapped edge of type 1, it will not have any other directly mapped edge. This 
follows from the definition of a directly mapped edge and Observation 12.11 There can be at 
most two directly mapped edges of the second tvpefObservation I2.2p . These edges mapped 
to e are always from a level < Level{e). There can be at most two directly mapped edges of 
type 3 also if they are not in T~L but are adjacent to e. By Observation 12.21 there can only be 
two such edges. 
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Figure 1: e G Ai. The edges marked * are not in % and are in M* . The edges which are not 
marked * are all in T~L. All the edges marked by * are indirectly mapped to e. 

Claim 2.2. There can be at most two directly mapped edges to an edge e £ Ai at any level. 

The total weight of the edges directly mapped to e will be maximum when both of them 
are from the same level as e. Assume that e is at level i. Summing the weights of the edges 
which are directly mapped to e, we get 

w(e*) < 2*a i+1 < 2aw(e) (1) 

Indirectly mapped edges can only be of the fourth kind in which the edge is not in H, 
but is adjacent to an edge in %, which in turn is adjacent to e. By definition, these edges are 
from a level lower than that of e. There can be at most two edges from each level lower than 
Level(e) which are in T~L and are adjacent to e(see Figure [1]). 

Claim 2.3. There can be at most two indirectly mapped edges to an edge e G Ai at level 
< Level(e). 

Note that there can be a large number of edges which are indirectly mapped to e. Still we 
will be able to get a good bound on their total weight. This is because there can be at most 
two indirectly mapped edges from each level and the weight of edges in the levels decreases 
geometrically as we go to lower levels. 

Assume that e is at level i. Summing the weight of edges which are indirectly mapped to 
e, we get 

Ei *\ o i+i 2a l+1 2aw(e) 
w e* < 2 ^ < < Y 2 

^-^ a — 1 a — 1 

e*G07 1 (e) J=*-l 

Thus, the total weight mapped to e is - 

£ w(e*)= £ w(e*)+ Y, w(e*)<w(e)(-^ rT + 2a\ 

e*e 9 i- 1 (e) e*e 9 i- 1 (e) e*e</»7 1 (e) 

As reasoned before, an edge in M* is mapped to some edge in Ai. So summing this over 
all the edges in At, we get 

eSM V ' eeM e *e</>- 1 (e) e*GM* 
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The function f(a) = y-^i + 2a J attains its minimum value of 8 at a = 2. So, if the value 
of a is picked to be 2, we get an 8 approximation maximum weight matching algorithm. We 
can state the following theorem. 

Theorem 2.1. There exists a fully dynamic algorithm that maintains 8-MWM for any graph 
on n vertices in expected amortized O(lognlogC) time per update. 

3 Improvements 

If e is an edge in M. at level i, then the approximation ratio calculated in the previous section 

can be written as 

_L y w(e *) = J_( V w(e*) + T w(e*) 
w(e) ^ w(e) V ^ ^ 



< 



e*e<£ z /(e) 6*60/(6) 

w(e*) \ 1 w(e*) 

2-^i 2 i ) 2i-3 ^ 2i 

Level(e*)=i Level(e*)=j 

For every edge e* at level i, if w{e*)/2 l < 7, then the approximation ratio is calculated as 



^ E -co s ( E t) + E^i E 7 

e*g</,- 1 (e) e'G^Ce) i =1 e*€0" 1 (e) 

Level(e*)=i Level(e*)=j 
i-1 

= 2 7 + V -^T.27 < 4 7 (By Claim E2] and E3D 

Z — / 2 l 3 

In the previous section, 7 = 2 in the worst case, so we get 8-approximation ratio. To 
improve the approximation ratio, we use randomization to find a good expected value of 7. 
To do this, we multiply all the weight of edges with a number chosen uniformly randomly 
from [1,2]. Since we are multiplying each edge with the same random number, M* does not 
change. Consider an edge e at level i having weight w{e) > 2 % before multiplication. Let X e 
be a random variable denoting the new weight of the edge e. Let Y e be a random variable 
denoting the relative weight of edge e. The relative weight of an edge e at level i is w{e)/2 l . 
Let w{e) = /32 m where 1/2 < j3 < 1. So X e lies in the range (f32 i+1 , f32 i+2 ). Y e is defined as: 

if X e e032*+i,2*+i], 
if X eG (2*+\/32*+ 2 ). 

Also, Y e is a continuos random variable having probability density function given by : 
f{x) = p2l +2 l _p T +i for all f32 i+1 < x < f32 i+2 The expected value of Y e gives the expected 
value of 7. The expected value of Y e is calculated as: 



Y 
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i / f 2i+1 X \ 1 / f 32 ^ 2 X \ 

E[Ye] = W+ 2 - P2 t+1 V V+i X; + P* +2 ~ /52 l+1 V J 2i+1 2^1 dX ) 

^ J 021+1 XAX + 2*+i y 2i+ i XdX J 



J TT (^r((2 J+1 ) 2 - /3 2 (2* +1 ) 2 ) + -^(/? 2 (2^ 2 ) 2 - (2^) 2 )) 



/32 
1 



I((l-/? 2 ) + (2/3 2 -l/2 



The function /(/?) = - (/3 2 + 1/2) has the maximum value when (3 = 1 and its maximum 

value is 1.5. From the above discussion, the expected value of 7 = 2£[Y e ]. This implies an 
approximation ratio of 4E[Y e ] = 6. 

Theorem 3.1. There exists a fully dynamic algorithm that maintains expected 6-MWM for 
any graph on n vertices in expected amortized O(lognlogC) time per update. 

Remark 3.1. For each 1 < r < 2, our algorithm computes a matching. This can be viewed 
as a mapping 4> from [1,2] to the set of matchings computed by our algorithm. The above 
analysis says that there exists a r* in [1,2] such that <p(r*) is an 6- approximation of MWM 
of the given graph. 

Note that our randomization step finds a number from [1,2]. We assumed that this random 
variable to be continuous over the above range. But in practice we only have O(logn) bits 
of randomness. In the appendix, we show that if we have clogn bits of randomness(c > 1), 
then we can get a approximation of (6 + e) where e = 0(l/n c ). 

4 Conclusion 

We presented a fully dynamic algorithms for maintaining matching of large size or weight 
in graphs. The algorithm for maintaining 6-MWM is the first fully dynamic algorithm for 
maintaining approximate maximum weight matching. 
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Appendix 

We show that if we have clogn random bits(c > 1), we get (6 + e) approximation of MWM. 
Assume that we get discrete random number between and 1 where the gap between two 
discrete points is l/n c . Consider the following function: f(i) = i/n c + 1 for < i < n c . The 
set {f(i)\0 < i < n c } gives the range of discrete random variables generated in the range 
[1,2]. Let e be an edge at level i having weight w(e) = f32 l+1 where 1/2 < f3 < 1. Let X r 
denote the new weight of an edge e. Given a value of r, the value of X r is f(r).(32 l+1 . When 
X r < 2 t+1 ,r < (L^J — l)n c and when X e = /32 l+2 ,r = n c . Also note that a discrete random 
number f(r) is generated with probability l/n c . Remember that Y e be a random variable 
denoting the relative weight of edge e. For the discrete case, we define it as follows: 

Ye \% ^ 0<r<( L ij-lK 

{^T, if (LjjJ -1K+1 <r<n c 

Expected value of Y e is calculated as: 



12 



E l Y e] = E jf-Pl /W generated ] + E 

r= ° ~ r=(^J-l)nc+l 



1 x — -\ .X^ 1 



(L^J-l)" 

r "° r= ( L lj_!) n c +1 

(L ^ l)nC /(r)/^ 1 1 ^ /(r)/?^ 1 1_ 

Z_> 2* 'n c ^ 2 i+1 ' n c 

r "° r=(L^J-l)nc+l 

(l£J-i)» c 







2 E /(*■) + E /(*■)) 

r= ° r=(L^J-l)n"+l 

= ^(E(^ c+1 )+ E (^ c +!)) 

< 4(S+- c + i )+( ( ^Vr n2c + ^-iK + i)) 

i?[Ye] has its maximum value when /3 = 1 and its maximum value is 1.5 + 2/n c . As 
discussed previously, the approximation ratio obtained by the algorithm is 4E[Y" e ]. This gives 
an approximation ratio of 6 + e where e = 8/n c . 
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